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[i] The amount of water stored and moving through the surface water bodies of large river 
basins (river, floodplains, wetlands) plays a major role in the global water and biochemical 
cycles and is a critical parameter for water resources management. However, the 
spatiotemporal variations of these freshwater reservoirs are still widely unknown at the 
global scale. Here, we propose a hypsographic curve approach to estimate surface 
freshwater storage variations over the Amazon basin combining surface water extent from a 
multi-satellite-technique with topographic data from the Global Digital Elevation Model 
(GDEM) from Advance Spacebome Thermal Emission and Reflection Radiometer 
(ASTER). Monthly surface water storage variations for ] 993-2007 are presented, showing 
a strong seasonal and interannual variability, and are evaluated against in situ river discharge 
and precipitation. The basin-scale mean annual amplitude of -1200 km 3 is in the range of 
previous estimates and contributes to about half of the Gravity Reco very And Climate 
Experiment (GRACE) total water storage variations. For the first time, we map the surface 
water volume anomaly during the extreme droughts of 1 997 (Qctober-November) and 2005 
(September-October) and found that during these dry events the water stored in the river 
and floodplains of the Amazon basin was, respectively, -230 (-40%) and 21 0 (-50%) km 3 
below the 1 993-2007 average. This new 1 5 year data set of surface water volume represents 
an unprecedented source of information for future hydrological or climate modeling of the 
Amazon. It is also a first step toward the development of such database at the global scale. 

Qlalioiu Papa, F., F. Frappart, A. Guntner, C. Prigent, F. Aires, A. C. V. Getirana, and R. Maurer (2013), Surface 
freshwater storage and variability in the Amazon basin from multi-satellite observations, 1993-2007, J. Geophys. Res. 
Atmos., 118 , 1 1,951-1 1,965, d oi: 1 0.1 002/201 3 JD020 500. 

1. Introduction 

[ 2 ] Terrestrial waters on Earth's icc-frec land represent less 
than 1 % of the total amount of water on Earth. However, they 
have a crucial impact on tores trial life and human environ- 
ment, and play a major role in climate vari ability. Excluding 
ice caps, fie sh water on land is stored in various reservoirs: 
snowpaek, glaciers, aquifers and other geological formations, 
root zone (upper few meters of the soil), and surface waters, 
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comprising of rivers, lakes, man-made reservoirs, wetlands, 
and inundated areas. These eomponents are continuously 
exchanging mass with the atmosphere and the ocean by verti- 
cal and horizontal motions, and thus are an integral part of the 
climate system with important links and feedbacks. Moreover, 
analysis of the flow, spatial distribution, and storage of 
freshwater on land are also a key issue for understanding the 
terrestrial branch of the global water cycle, as well as for the 
management of water resources [Chahine, 1 992; Kundzewicz 
et al . , 2007]. 

[ 3 ] However, although improved description of the com- 
ponents of the global water cycle is now recognized as being 
of major importance, the global distribution and spatiotem- 
poral variations of continental water extent and volume arc 
still poorly known [Alsdorf et al ., 2007a; Papa el al., 20 1 Oaf 

[ 4 ] Until recently, our knowledge of the spatiotemporal 
variations of continental waters relied on sparse in situ obser- 
vations and hydrological models. In situ gauge measurements 
help quantify the movement of water discharge in river chan- 
nels but, for instance, provide comparatively little information 
about the spatial dynamics of terrestrial water in floodplains 
and wetlands or groundwater. In addition, the availability of 
ground-based gauge information has dramatically decreased 
during the last decades [Alsdorf and Lettenmaier , , 2003; 
GRDC, Global Runoff Data Base-Statistics 2012, http:// 
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www.bafg.de/GEU3C/EN/01 GRDC/13 dtbsc/db stat.html? 
nn-201876], especially in remote areas with difficult access 
such as tropical regions. In some eascs 5 public access to the 
latest river discharge observations can also be restricted. 
Lacking spatially complete measurements of inundation/wct- 
land locations, sizes, and water volume changes, hydrologic 
models arc unable to properly partition precipitation (minus 
evapo transpiration) among these several components and 
represent their effects on river discharge at contincntal-to- 
global scales [Alsdorf and Lettenmaier, 2003; Alsdorf et al., 
2007a; Decharme et al., 2008, 2012; YamazaM et al., 2011; 
Getirana et al., 2012]. Our ability to measure, monitor, and 
forecast supplies of freshwater using in situ methods and 
hydrological models is facing considerable difficulties, at least 
at large scales. 

[s] The terrestrial water balance equation [Peixoto and 
Oort , 1992], as applied to a river basin, can be written as 

P-E = QsPQg + AV, (1) 

with P and E representing basin-averaged precipitation and 
evapotranspiration, respectively; Qs is river discharge; Qg 
is groundwater discharge across the basin boundary; and V 
is the total surface and subsurface storage (sum of soil mois- 
ture, snow water content, surface water storage, vegetation wa- 
ter content, groundwater, and glaciers [. Lettenmaier , 2005]), In 
most eases, each quantity is available only with a large uncer- 
tainty [Sheffield et al., 2009, Azarderakhsh et al., 2011] and 
the estimation of the storage change term is particularly prob- 
lematic. Given that -60% of the world floodplains and wet- 
lands arc only inundated at some time in the year [Matthews, 
2000], intraseasonal and interannual variations in surface- 
stored water volumes at regional -to-global scales arc not well 
known and limit our ability to understand their impacts on 
evapotranspiration, infiltration, and runoff. With the absence 
of a coherent observational basis for quantifying water stor- 
age globally over long time periods, the water storage term 
is often ignored. However, this approach precludes consider- 
ation of surface and subsurface water volume dynamics at 
shorter time scales. 

[ft] Remote sensing techniques have been very useful to hy- 
drology investigations especially over the last 20 years [Smith, 
1997; Cazenave et al., 2004; Alsdorf and Lettenmaier, 2003, 
Alsdorf et al. , 2007b; Hess et al. , 2003; Prigent et al ., 2007; 
Papa et al., 2010b; Papa et al., 2012a, among others]. 

[7] The concept of measuring the hydraulics of inland 
water bodies from space was first brought forth in the late 
1 990s based on the suceesses ofthe Topcx/POS EICON radar 
altimetry mission that provided the systematic monitoring of 
water levels of large rivers, lakes, and floodplains [Birkett 
et al., 2002; Cretaux et al., 2005]. Multi-satellite remote 
sensing techniques [e.g., Prigent et al., 2001, 2007, 2012; 
Papa et al ., 2006, 2008a, 2010a] also offer important infor- 
mation on land surface waters, such as the variations of 
surface water extent at the global scale. Since 2002, the 
Gravity Recovery And Climate Experiment (GRACE) grav- 
ity mission offers, for the first time, direct estimates of the 
spatiotcmporal variations of Total terrestrial Water Storage 
or TWS (the sum of ground water, soil water, surface water, 
and snowpack) [Ramillien et al., 2005; Rodell et al., 2007; 
Tapley etal., 2004] at seasonal and basin scales. Nevertheless, 
despite the increasing number of satellite observations, we still 


have surprisingly poor knowledge of the dynamic of surface 
freshwater storage at least globally and at large scale. Lacking 
such observations at proper space and time scales, several ba- 
sic questions related to the land surface water budget are still 
open; perhaps the most fundamental being: flow much fresh- 
water is stored at the surfkcc/ncar surface of continents? 
What are the spatial and temporal dynamics in terrestrial sur- 
face water storage and their link with climate variability? 

[s] Recently, some efforts have been undertaken to quan- 
tify the surface freshwater storage and its variations at sea- 
sonal to interannual time scale using satellite observations. 
The technique, which combines surface water extent obser- 
vations from a multi-satellite technique [Papa et al., 2010a; 
Prigent et al., 2007, 2012; the Global inundation Extent from 
Multi-Satellite called hereafter (CiEMS)] and radar altime- 
ter-derived height variations of rivers, wetlands, and inunda- 
tions [Frappart etal., 2006], was first developed over the Rio 
Negro, a subbasin of the Amazon [Frappart et al., 2008, 
2011a], and also tested over a boreal environment in the 
Ob River basin [Frappart et al., 2010a]. Using continuous 
water level observations derived from ENVISAT radar altim- 
eter [Sa/xfos da Silva et al., 2012] between 2003 and 2007, 
Frappart et al. [2012] provide for the first time monthly var- 
iations of surface water storage for the entire Amazon basin. 
In particular, the results highlighted that during the 2005 ex- 
ceptional drought, the amount of water stored in the river and 
the extensive floodplains (covering more than 300,000 km 2 
or 5% of the surface of the entire basin) was —130km 3 
(-70%) below its 2003 2007 average. 

[ 9 ] The combination of GIEMS and altimeter observations 
has proved to be a new powerful and reliable tool for moni- 
toring large-scale surface freshwater dynamics. However, 
the availability of a large sample of altimctry-dcrivcd water 
levels (more than 530 virtual stations for the Amazon basin 
case study) used in Frappart et al. [20 1 2] remains an excep- 
tion. Indeed, former and current radar altimeter missions, char- 
acterized by an orbit intertrack interval of- 80 km (ERS-1/2, 
ENVISAT, AltiKa) and -300 km (T/P and Jason-2) at the 
equator, do not provide an adequate monitoring with sufficient 
spacc/timc sampling of most of continental water bodies 
worldwide [Bianeamaria et al., 2010; Papa et at , 2012b] 
and thus limit, at the time of writing, the development of this 
technique at global scale. 

[iu] In the present paper, we propose an innovative 
observation -based technique to estimate surface freshwater 
storage variations in rivers, floodplains, lakes, and wetlands. 
It combines the surface water extent from GIEMS with topo- 
graphic data derived from the Global Digital Elevation 
Model (GDEM) from the Advance Spacebomc Thermal 
Emission and Reflection Radiometer (ASTER), using a 
hypsographic curve approach. Keeping in mind that the ulti- 
mate objective is to derive the spatial and temporal variability 
of freshwater storage at the global scale, the present study 
proposes to focus on the development of the new technique 
and on the analysis of the results over the Amazon River ba- 
sin, the largest hydrological system of the world (-6 million 
km 2 of surface area), characterized by extensive floodplains 
and contributing with 15-20% to the total freshwater dis- 
charge to the world oceans [Richey et al., 1 986]. 

[ 11 ] Section 2 will present the data sets used in this study. 
In section 3, focusing on the Amazon basin, we will present 
the methodology of the hypsographic curve approach, and 
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the combination of G1EMS with ASTER GDEM data set 
to derive surface water storage. In section 4, the results arc 
presented and discussed over the 1993-2007 period. An 
evaluation is performed comparing the new estimates with 
other external data sets such as GRACE-dcrivcd TWS, in situ 
river discharge observations, and precipitation. Analysis of 
the results is also presented with emphasis on the two excep- 
tional droughts of 1997 and 2005. Finally, conclusions and 
perspectives are discussed in section 5. 

2_ Data Set 

2-1. Multi-Satellite-Based Surface Water Extent Dynamic 

[ 12 ] The complete methodology which captures at the 
global scale the extent of episodic and seasonal inundations, 
wetlands, rivers, lakes, and irrigated agriculture is described 
in detail in Prigent et al. [2001, 2007, 2012] and Papa 
et al. [2006, 2010a]. The technique uses a complementary 
suite of satellite observations covering a large wavelength 
range: (1) Advanced Very High Resolution Radiometer 
(AV11RR) visible (0.58 0.68 pm) and near-infrared (0.73- 
1.1 pm) reflectances and the derived Normalized Difference 
Vegetation Index (NDVX), (2) passive microwave emissivi- 
ties between 19 and 85 GHz. They are estimated from the 
Special Sensor Microwave/Imager (SSM/1) observations 
by removing the contributions of the atmosphere (water 
vapor, clouds, rain) and the modulation by the surface tem- 
perature [Prigent et al. , 1997; Prigent et al. , 2006]. The tech- 
nique uses ancillary data from the International Satellite Cloud 
Climatology Project (ISCCP) [Roraow and Schiffer, 1 999] and 
the National Centers for Environment Prediction (NCEP) 
rcanalysis [Kalnay et al., 1996], and (3) backscattcr at 
5.25 GHz from the European Remote Sensing (ERS) satel- 
lite scattcromctcr. 

[n] Observations arc averaged over each month and 
mapped to an equal area grid of 0.25° resolution at the equa- 
tor (each pixel equals 773 km 2 ) [Prigent et al. , 200 1 ; Prigent 
et al ., 2006]. An un supervised classification of the three 
sources of satellite data is performed and the pixels with 
satellite signatures likely related to inundation arc retained. 
For each inundated pixel, the monthly fractional coverage 
by open water is obtained using the passive microwave signal 
and a linear mixture model with end-members calibrated with 
scattcromctcr observations to account for the effects of vege- 
tation cover [Prigent et al. , 2001, 2007]. As the microwave 
measurements arc also sensitive to the snow cover, snow 
and ice masks arc used to filter the results and avoid any con- 
fusion with snow-covered pixels [Armstrong and Brodzik, 
2005]. Because the ERS scatteromcter encountered serious 
technical problems after 2000, the processing scheme had 
to be adapted to extend the data set and monthly mean clima- 
tology of ERS and NDVl-AViiRR observations arc used 
[Papa et al. , 2010a; Prigent et al ., 2012]. Fifteen years of 
global monthly water surfaces extent 1 993-2007 is available 
[Prigent et al ., 2012]. The data set has been extensively eval- 
uated at the global scale [Papa et al., 2010a; Prigent et al., 
2007; Papa et al ., 2008a] and for a wide range of environ- 
ments, including boreal regions [Papa et al., 2007, 2008b] 
and tropical regions [Papa et al., 2006; Frappart et al., 
2008]. This data set is intensively used for climatic and hy- 
drological analyses, such as the evaluation of methane sur- 
face emissions models [Bousquet et al., 2006; Ringeval 


et al., 2010; Wania et al., 2013] and the validation of the 
river flooding schemes coupled with land surface models 
[Decharme et al., 2008, 2012; Ringeval et al., 2012; Getirana 
et al., 20 1 2; Pedinotti et al., 2012]. 

[ 14 ] Figure 1 shows GIEMS characteristics over the Amazon 
River basin (Figure la). Figures lb and le show, respectively, 
the annual mean and annual maximum extent of surface water 
averaged over 1 5 years (1 80 months). They exhibit very realis- 
tic distributions of major rivers and tributaries (Amazonas, 
Solimoes, Negro) with associated inundated areas and wet- 
lands well delineated even in complex areas characterized by 
extensive flooding below dense vegetation canopies, such 
as the Central Amazonian corridor. The spatial distribution of 
G1EMS over the Amazon basin was extensively evaluated 
against high-resolution (100 m) SAR images [Hess et al., 
2003] in Prigent et al. [2007] and in Aires et al. [2013], as well 
as using other regional surveys representing various compo- 
nents of wetland and open- water distributions. The evalua- 
tion led to an overall estimation of GIEMS uncertainties of 
-1 0%. Seasonal and interannual variations of the basin-scale 
total surface water extent and the associated anomalies are 
presented in Figures 1 d and 1 e. The extent has a strong sea- 
sonal cycle (Figure Id), with a mean annual averaged maxi- 
mum of~2.5 x 10 5 km 2 for 1 993 2007. It shows a substantial 
interannual variability especially near the minima and max- 
ima. The years 1997, 2002, 2006, and 2007 exhibit larger 
peaks (Figures Id and le) mostly associated with El Niiio- 
La Nina events. During El Nino years, the eastern part of 
the Amazon basin is drier than normal years, whereas the 
Andean catchment at medium and high elevation presents a 
slight increase in rainfall. During la Nina years, the south- 
western parts of the Amazon basin receive more rainfall 
than normal years, especially in the Altiplano and high eleva- 
tion regions of the Eastern Cordillera [Espinoza et al., 2009; 
Bookhagen and Strecker, 2010]. The droughts of 1997 and 
2005 [Marengo et al., 2008a; Tomas el la et al., 2011; 
Frappart et ai, 2012] arc also well captured. 

2.2. ASTER GDEM 

[ 15 ] ASTER GDEM was developed jointly by the National 
Aeronautic and Space Administration (NASA) and the 
Ministry of Economy, Trade, and Industry (MET1) of Japan. 
The ASTER instrument, launched onboard NASA's Terra 
spacecraft in December 1999, has along-traek stereoscopic 
capability using its near infrared spectral band and its nadir- 
viewing and backward- viewing telescopes to acquire stereo 
image data with a base-to-height ratio of 0.6. The basic charac- 
teristics of stereoscopy and its application to the ASTER sys- 
tem for GDEM generation are explained in detail in Toutin 
[2008]. The horizontal spatial resolution is 1 5 m and one na- 
dir-looking ASTER visible and near-infrared (VN1R) scene 
corresponds to about 60 x 60 km ground area. The methodol- 
ogy used to produce the ASTER GDEM involves automated 
processing of the entire 1 .5-million-sccne ASTER archive 
acquired from the start of observation until August 2008 
[Toutin, 2008; Abrams et al., 2010]. The processing includes 
stereo-correlation to produce -1,264,000 individual scene- 
based AS TER DEMs, cloud masking to remove cloudy pixels, 
stacking all cloud-screened DEMs, removing residual bad 
values and outliers, averaging selected data to create final pixel 
values, and then correcting residual anomalies. The ASTER 
GDEM covers land surfaces between 83°N and 83°S and is 
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Figure I. The global inundation extent from multi -satellite (GIEMS) and anomalies over the Amazon basin, 
(a) The Amazon basin and its major rivers and tributaries. The red star shows the locations of the gauging 
station ObidoSj Para, Brazil, located 8001cm upstream from the mouth of the Amazon, (b) and (c) Map of 
annual mean and annual maximum surface water extent averaged over 1993- 2007, for each 773 km 2 pixel 
(d) Monthly mean surface water extents for 1 993- 2007 for the entire Amazon basin, (e) Corresponding anom- 
alies (by removing the mean 1993-2007 (red line) and deseason alized obtained by subtracting the 1 5 year 
mean monthly value from individual months (black line)). 


partitioned into 22,600 tiles of 1° by 1° (containing at least 
0.01% of land area). ASTER GDEM has a 1 arc-second 
(30 m) grid of elevation and is referenced with respect to the 
WGS84/ECM96 gcoid. Several studies have dealt with the 
evaluation of ASTER GDEM at local to regional scales 
[Hirano et aL, 2003; Fujisada et aL, 2005; Hayakawa et aL, 
2008; Peng et aL, 2013]. Preproduction accuracies for the 
global product were estimated at --20 m at 95% confidence 
level for vertical data and -30 m at 95% confidence level for 
horizontal data. In this study, we use the ASTER GDEM 
Version 2 released in October 201 1 available at http://gdem. 
ersdae .j spacesy stems, or. j p/search. j sp . 

2.3. Complementary Data for Evaluation 
2.3.1. GRACE Data 

[i6] The Gravity Recovery And Climate Experiment mis- 
sion, launched in March 2002, provides measurements of 
the spatiotemporal changes in Earth's gravity field. Several 
recent studies have shown that GRACE data over the conti- 
nents can be used to derive the monthly changes of the total 
land water storage [Ramil lien et aL, 2005; Schmidt et aL, 
2008; Lander er and Swenson, 2012] with an accuracy of 
~1.5 cm of equivalent water thickness when averaged 
over surfaces of a few hundred square-kilometres. In this 
study, we use monthly solutions from the Geo Forsehung 
Zentrum (GFZ), the University of Texas Center for Space 
Research (UTCSR), and the Jet Propulsion Laboratory 
(JPL) from February 2003 (data are missing for January 
2003) to December 2007 in order to analyze the time varia- 
tions of the water mass changes in the Amazon basin. 
Unfortunately, the GRACE solutions suffer from the 


presence of an unrealistic high-frequency noise correspond- 
ing to north-south striping that is caused by orbit resonance 
during the Stokes coefficients determination and aliasing of 
poorly modelled short-term phenomena. To attenuate the 
noise in the Lcvel-2 GRACE solutions, we used the global 
solutions postprocessing by an Independent Component 
Analysis (1C A) approach based on the combination of GFZ / 
UTCSR/JPL solutions of the same monthly period to isolate 
statistically independent components of the observed grav- 
ity field, and mainly the continental water storage contribu- 
tion [Frappart et aL, 2010b, 2011b]. These data can be 
downloaded at: http://grgs.obs-mip.fr. 

2.3.2. River Discharge 

[ 17 ] Monthly discharges observed at Obidos (Figure la), 
state of Para, Brazil, which is the closest gauge to the mouth 
of the Amazon River (-800 km), is used in this study. The re- 
cord is available at the Environmental Research Observatory 
(ORE) Gcodynamieal, hydrological, and bfogcochcmical 
control of crosion/alteration and material transport in the 
Amazon basin (HYBAM) website (http ://www.orc-hy bam, 
org/) over the period 1993-2007. 

3, Methodology 

3.1. The Ilypsographic Curve Approach 
With ASTER GDEM 

[is] The method to estimate surface freshwater storage 
consists of the combination of the surface water extent from 
GIEMS product with ASTER GDEM, using a hypsographic 
curve approach that relates the flooded area to the elevation. 
It is a three-step process (Figure 2) that can be summarized 
as follows: 
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Figu re 2, Schematic representation of the hypsographic curve 
technique to estimate surface water volume horn GIEMS (see 
text for details). 

[19J I. For each cell of the GIEMS data set, the cumulative 
distribution function of elevation values is first derived from 
the corresponding subset of ASTER GDEM. For each pixel 
of the surface water extent data set (equal -area grid of 
773 km 2 ), the corresponding subset of the ASTER GDEM 
is selected. Figure 3a, 3c, 3e, and 3g (left panel) shows 
the histograms of ASTER GDEM elevation within tour 
pixels (equal-area grid of 773 km 2 ) at different locations 
in the Amazon basin, located in the Marnore floodplain 
(Figure 3a), the central Amazonian floodplain (Figure 3c), 
the Tapajos floodplain (Figure 3e), and the Rio Negro flood- 
plain (Figure 3g). Approximately 800,000 elevation points 
of ASTER DEM fall within the satellite-derived surface wa- 
ter extent cell. The next step is to construct the so-called 
hypsographic curve, or curve of cumulative frequencies, for 
each cell of GIEMS data set based on ASTER GDEM eleva- 
tions. The hypsographic curve is then equivalents the distri- 
bution of elevation values in each 773 km 2 cell sorted in 
ascending order and represents an area-elevation relationship 
as illustrated in Figures 3b, 3d, 3f, and 3g (right panel) corre- 
sponding to the four locations in the Amazon basin. 

[ 20 j 2. First, for each pixel of GIEMS, a translation is 
applied to set to zero the lowest elevation of the hypsographic 
curve by subtracting the lowest value from all other eleva- 
tions. The hypsographic curve is then converted into an 
area-surface water volume relationship by estimating the sur- 
face water volume associated w ith an increase of the pixel 
fractional open water coverage by filling the hypsographic 
curve from its base level to an upward level, following; 

V(U) = !(*(/) h{i 1)) x S/100 x i (2) 

i—l 

where V is the surface water volume (in km 3 ) for a percentage 
of flood/in undation □ (a step of 1 % in increment i in percent- 
age of inundation is chosen here), S the 773 km 2 area of a 
GIEMS pixel, and h the elevation (in km) for a percentage 
of flood/inundation C given by the hypsographic curve. 
Examples of surface-volume profile are shown in Figure 4 
for 16 ioeations around the Amazon basin (wc will call 
without distinction “hypsographic curve,” the area-elevation 


relationship as in Figure 3 or the area-volume relationship as 
in Figure 4). 

[ 21 ] 3. Finally, in the last step, in order to estimate the sur- 
face water volume, the hypsographic curve (as in Figure 4) 
obtained tor each grid cell of the Amazon basin is combined 
with the monthly variations of surface water extent from 
GIEMS. The surface water volume tor each pixel and each 
month is estimated by intersecting the pixel hypsographic 
curve with the GIEMS estimates of pixel water coverage 
for that month, as shown in Figure 2. Note that, with the pro- 
posed method, the water storages below the lowest levels of 
storage are not accessible. Thus, the estimated water storage 
represents the increment above the minimum storage. 

[ 22 ] Figure 4 presents, for i 6 locations in the Amazon basin, 
the comparison between the ASTER GDEM hypsographic 
curves and the ones constructed by Frappart et al. [2012J. 
The hypsographic curves from Frappart et al. [2012] were 
obtained using GIEMS surface water bodies in combination 
with water level variations derived from ENVISAT Radar 
Altimeter RA-2 (processed using the lce-1 retracker as in 
Santos da Silva et al . [2012]). Note that botli hypsographic 
curve data sets have been derived in a totally independent 
way. Figure 4 clearly shows, for most cases, a very good 
agreement between both hypsographic curve approaches. 



Elevation (m) Cumulative Area (% al pixel surface) 


Figure 3. Examples of hypsographic curves from ASTER 
GDEM for the Amazon basin. Left: The histogram of ASTER 
GDEM elevation for a 0.25° equal grid area. Right: The 
hypsographic curves from ASTER GDEM, i.e., the distri- 
bution of elevation values in each 773 km 2 cell sorted in 
ascending order. 
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Figure 4. Surface volume profile (equivalent to the hypsographie curve), i.e., the relationship between 
surface water storage within each grid cel! and the inundated area of a 773 km 2 pixel (in percent). Red 
curves are derived from ASTER GDEM and black curves are obtained using ENVTSAT radar altimeter 
as in Frappart et al. [2012]. The blue (green) dashed line is the maximum (minimum) coverage of surface 
water observed by GIEMS during 1993-2007. 


Small differences in water volume, in general less than 
Q.5km 3 , are observed for surface water extent covering 
up to 60-70% of a pixel (Figures 4a, 4c, 4d, 4c, 4j, and 
4m) and up to 80 90% of a pixel (Figures 4g, 4i, 4k, 41, 
and 4o). Some eases show a good agreement in the low wa- 
ter coverage range (<40% coverage) but a mismatch for 
higher flood area extent (Figures 4b, 4f, and 4h). Some cases 
show a larger mismatch such as in Figures 4n and 4p. 
However, one can notice a general behavior with ASTER 
GDEM hypsographie curves: potentially large fractions of 
surface water extent correspond to unreasonably high surface 
water storage (note that to make the figures comparable, the 
range of water volume values was limited to 20 km 3 , but 
some values at the higher end of the curve can reach up to 
100 km 3 and higher). Some eases (Figures 4a, 4c, 4e, 4g, 4i, 
4j, 4k, 41, 4m, 4n, and 4o) even show an extremely sharp in- 
crease in surface water volume for the last 20% coverage of 
the pixel, going from -5 to 10 km 3 for 80 90% coverage to 
hundreds of km 3 for 100% of flood. This is not observed in 


the ENViSAT-dcrivcd hypsographie curves for which in 
general the maximum in surface water volume is potentially 
-10 to 15 km 3 . The blue and green dashed lines in Figure 4 
give, respectively, the maximum and minimum surface water 
extent observed by GLEMS during the period 1993^2007. It 
shows that in most eases the sharp increase of the surface vol- 
ume profile around the higher end points will not affect the 
estimates as such high flood extents are never reached. 
Nevertheless, when a large surface water extent at maximum 
is observed with GIEMS, such as in Figures 4c and 4k and 
especially Figure 4n where 100% of flooding is found in 
the main channel of the Amazon, it can lead to large overes- 
timations of surface water volume. In these eases, values are 
jumping from -6 kin 3 at 90% to -10 km 3 at 95% and to 
an unrealistic -230 km 3 at 100% coverage of flooding for a 
single pixel. As a result, applying the step 3 of the methodol- 
ogy for all pixels of the Amazon basin leads to an unrealistic 
surface water volume with a mean annual amplitude of 
-3500 km 3 . If represents more than 1.5 times the seasonal 
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Figure 5. For the same locations as in Figure 4, the standard deviation (STD, in km 3 ) of the water volume 
derived from ASTER GDEM hypsographic curve (red plus sign) calculated over 5% flood coverage win- 
dows (20 windows named Nb on the x axis) calculated for each curve, see text for details). Same for 
ENVISAT-derived hypsographic curve with black plus sign. 


amplitude measured by GRACE (-2300 km 3 ), composed of 
surface water volume, soil moisture, and root zone and ground- 
water storage, In FrappartetaL [2012], the seasonal amplitude 
of the Amazon surface water storage was found to be of 
—1000 km 3 , which lies within the same range as large-scale 
model simulations. For instance, Get ir ana etai [201 2] reported 
a seasonal amplitude of -1100 km 3 using the Hydrological 
Modeling and Analysis Platform (IlyMAP) model, which 
includes floodplain reservoirs, whereas surface water storage 
variations were found to be of -470 km 3 when simulated by 
the WaterGAP Global 1 lydrology Model (WG11M) with a very 
simple floodplain storage module [Giintner et ai , 2007]. 

[ 23 ] The hypsographic curves showing an extremely sharp 
increase in surface water volume result from ASTER GDEM 
uncertainties and potential errors. Indeed, one of the major 
limitations of satellite-derived OEMs is that they arc not al- 
ways representing bare earth but can include vegetation and 
man-made structures. Moreover, because the acquisition of 
elevation data with ASTER is made using near infrared spec- 
tral band, GDEM can be affected by cloud cover, such as 
very low but dense boundary layer clouds in tropical regions. 


These data arc difficult to filter in the raw ASTER GDEM de- 
spite large processing and can result in erroneous high 
elevation topographic data, inducing further large errors in 
the hypsographic curves. 

[24] In order to avoid spurious contamination of ASTER 
GDEM-based surface water volumes due to high elevation 
values of some pixels at the upper edge of the hypsographic 
curve, we propose in the following a general correction of 
ASTER GDEM hypsographic curves. 

3.2, Correction of ASTER GDEM Hypsographic Curves 

[ 25 ] Figure 5 shows, for the same locations as in Figure 4, 
the Standard Deviation (STD) of the water volume derived 
from ASTER GDEM hypsographic curve (red ‘‘plus” sign) 
calculated over 5% flood coverage windows (20 STDs calcu- 
lated for each single curve). For instance, the first plus sign 
represents the STD of the surface water volume calculated 
using the values at 1 , 2, 3, 4, and 5% (first window, Nb = 1 
in Figure 5), the second plus sign represents the STD calcu- 
lated with the corresponding values at 6, 7, 8, 9, and 10% 
(second window, Nb = 2), and so on. The STD values 
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Figure 6. Examples of surface volume profile (similar as in Figure 4), i.e., the relationship between sur- 
face water storage within each grid cell and the inundated area for the Amazon basin. Red curves are de- 
rived from ASTER GDEM, black curves are obtained using ENVISAT radar altimeter as in Frappari 
et aL [2012], and green curves are derived from ASTER GDEM after correction (red curves corrected). 


estimated from ENVfSAT-dcrivcd hypsographie curves are 
also shown with a black plus sign. Figure 5 shows that in 
most cases, except for few isolated values, the STD calcu- 
lated with ENVISAT hypsographie curves utc below 
0.3 km 3 . it means, in other words, that within 5% of a 
773 km 2 pixel, i.e., -35 km 2 , the surface water volume does 
not vary more than 0.3 km 3 , corresponding nevertheless to 
water level variations of— 8.5 m at maximum. These are real- 
istic values of volume and water level variations observed in 
the Amazon main channel and floodplains [F rap part et a /. , 
2006; Santos Da Silva et aL, 2012; Getirana et aL , 2012]. As 
expected from Figure 4, STD values calculated from ASTER 
GDEM hypsographie curves match well with ENVISAT 
values, especially for the 1 5 first STD values corresponding 
to —75% flood coverage. However, Figure 4 also shows that 
for most cases, the STD values calculated from ASTER 
GDEM hypsographie curves for coverage above 80% are 
much higher than the ENVISAT ones. For ASTER GDEM, 
several consecutive values of STD are found between 0.6 and 
1 km 3 and above. These values are mostly unrealistic as a 


variation of 0.8 km 3 in water volume corresponds to a water 
level increase of more than 22 m by 35 kin 2 increment. 

[ 26 ] We propose a simple procedure to correct the behav- 
ior of ASTER GDEM hypsographie curves to prevent 
overestimation of surface water volume at pixel level. 
For each percent value of flood coverage area (see 
Figure 4), if the value belongs to a 5% window where the 
STD is below a 0.4km 3 threshold, then the corresponding 
surface water volume is kept. Subsequently, if the percent 
value belongs to a 5% window where the STD is above 
0.4 km 3 , the corresponding surface water volume is re- 
placed by the fitted value based on a simple linear regres- 
sion analysis using the 10 previous water volume values 
of the hypsographie curve. For example, a given point at 
80% flood coverage, which belongs to a window with a 
STD more than 0.4 km 3 , will be replaced by the fitted value 
computed using the simple linear regression equation 
obtained from the values between 70 and 79%. 

[ 2 ?] Figure 6 shows the results, for the same locations as in 
Figure 4, after the correction is applied to all AST ER GDEM 
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Figure 7. Surface water volume of the Amazon, (a) Monthly mean surface freshwater volume variations 
for 1993- 2007 for the entire Amazon basin as derived from the combination of G1EM3 estimates and 
ASTER GDEM hypsographic curves (black curves) and compared to Total Water Storage variations esti- 
mated from GRACE (green), (b) Mean seasonal cycle (2003-2007) of Amazon surface freshwater volume 
variations (black) and GRACE Total Water Storage (green). 


hypsographic curves (green line). It shows in most eases 
more realistic values even for large flood extent. For each 
of the 16 locations shown in Figure 6, the Root Mean Square 
Difference (RMSD) between EN VISAT-derived hypsographic 
curves and ASTER-derived original hypsographic curves 
ranges from 3.39 km 3 (Figure 6i) to 33.8 km 3 (Figure 6p). 
After correction, the RJV1SD between EN VISAT-derived 
hypsographic curves and ASTER-derived hypsographic curves 
is reduced and ranges between 0.34 km 3 and 1 .05 km, respec- 
tively. For all 16 locations, the mean RMSD between 
EN VISAT-derived hypsographic curves and ASTER-derived 
original hypsographic curves is —10.1 knr whereas it is 
-0.62 km 3 when using ASTER-derived hypsographic curves 
after correction. 

[ 28 ] When looking at the low end of the hypsographic 
curves, it should be also noted that the proposed method 
assumes that we cannot have access to water storage below 
the minimum values that ASTER GDEM (and the G1EMS) 


can provide. This can be a potential source of error when 
estimating the extreme low storage values of exceptional 
drought years. Indeed, in order to capture correctly the ex- 
treme low storage values during droughts, ASTER GDEM 
should have produced credible elevation data for those pe- 
riods at the low end of the histograms. Unfortunately, it is 
not possible to verify such information. 

4. Results, Evaluation, and Discussion 

[ 29 ] Combining the corrected AST ER GDEM hypsographic 
curves and the G1EMS satellite-derived observations (Step 3), 
wc can now estimate for the first time the long-term 
Amazon surface water storage and variations for the period 
1993 2007. Figure 7 shows the monthly variations 1993 
2007 (Figure 7a, 1 5 years) and seasonal cycle (calculated 
over 2003-2007, Figure 7b) of the surface water volume 
aggregated for the entire Amazon basin and compared to 
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Figure 8. Maps of Amazon surface freshwater volume from 1 993 to 2007. (a) Annual mean surface water 
volume averaged over 1993 2007, for each 773 km z pixel, (b) Annual maximum surface water volume 
averaged over 1993 - 2007. (c) The difference between the annual maximum and minimum surface water 
volume averaged over 1993 2007. (d), (e) and (f) Same as Figures 8a, 8b, and 8c but zoomed over the 
Amazon main corridor. 
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Figure 9* Correspondence among the satellite-derived monthly mean surface water storage, in situ river 
discharge and GR ACE-based Total Water Storage in the Amazon, (a) For the entire Amazon basin, the 
black line is the satellite-derived surface water storage, the green line is the GRACE-based Total Water 
Storage, the red line is the in situ river discharge from the Brazilian Water National observed at Obidos. 
(b) Mean seasonal cycle (2003-2007) of Amazon surface (black) and total freshwater storage variations 
(green) and in situ river (red)- (c) Deseasonalized normalized anomalies (obtained by subtracting the 
1 5 year mean monthly value from individual months and dividing by the standard deviations of the raw 
time series) for satellite-based surface wafer storage (black line) and in situ river discharge (red line). 
The deseasonalized normalized anomalies of basin-scale mean precipitation from GPCP is also shown 
(blue line), (d) Corresponding deseasonalized anomalies tor 2003-2007 along with the GR ACE-derived 
total water volume change anomalies (green). The purple solid line is the deseasonalized anomalies for 
2003-2007 of the surface water volume (multiply by 5 tor visual purpose) aggregated over a large region 
[0°S-4°S; 57°W— 60°W] west of Obidos location, 


GRACE-dcrived TWS. In parallel, Figure 8 shows the spa- 
tial distribution (annual mean, annual mean maximum, and 
mean annual amplitude) of surface water volume for the 
entire basin with a focus on the Amazon central corridor. 

[ 30 ] Results show a strong seasonal cycle (Figures 7a and 
7b), with a mean annual amplitude of ~1 200 km 3 with a max- 
imum surface water volume observed in March, two months 
ahead of GRACE TWS annual peak, and a minimum value 
reached in November, one month after GRACE TWS mini- 
mum is recorded. This delay can be explained by the slower 
groundwater flow in comparison to the surface water move- 
ment, causing the 2 month delayed peak in the GRACE 
TWS time series. Figure 7b also shows that the surface water 
reservoir variation represents about 50% of the TWS varia- 
tions as measured by GRACE (Figure 7b). This is in the same 
order of magnitude as previous accepted results on the parti- 
tion of TWS into contributing hydrological storages [Pokhrel 
et al. , 2013]. For instance, using remote-sensing observa- 
tions over the Rio Negro, one of the major tributaries of 
the Amazon River, [Frappart et al. , 2008; Frappart et al. 
2011a] showed that the surface water storage contributes to 
50 60% of the TWS variations. At basin scale, only results 
based on model simulations arc available. For instance, using 
hydrological model outputs based on a framework combin- 
ing a land surface model and a global runoff routing scheme, 
Kim et al. [2009] suggested that river storage explains 
~73% of TWS variation. Han et al. [2009], based on simu- 
lations from Global Land Data Assimilation System (GLDAS) 
Noah Land Surface Model and Alkama et al. [2010], using 
the interactions between Soil-Biosphere- Atmosphere (1SBA) 
land surface model, both indicated that TWS variations in 
the Amazon are almost equally partitioned into soil moisture 
and river storage variations. Modeling results from WG11M 
[Gimtner et al., 2007] suggested that surface water storage 
contributes to -40% of seasonal TWS variations in the 


Amazon. Recently, Paiva et al. [2013] showed that using the 
new MGB-iPH model, surface waters dominate TWS for the 
whole Amazon area with a fraction of 56%. 

[si] The new estimates of surface water storage also show 
a substantial interannual variability at basin scale, especially 
in terms of annual maximum and minimum. For instance, the 
years 1997, 2002, 2006, and 2007 exhibit larger peaks and 
the years 1 997 and 2005 show more extreme minima associ- 
ated with the major drought during these years (discussed 
later). The mtcrannual variability for the years 2003-2007 
(60 months) is also in close agreement with the inter annual 
variations observed in the GRACE-dcrived TWS, with a cor- 
relation coefficient /? _ 0.95 (p-valuc < 0.01) and a lag-time 
of 2 months. 

[n] In terms of spatial distribution, Figure 8 shows realis- 
tic structures with larger surface water volume and changes 
observed in the main Amazon channels and major tributaries 
(Negro, Solimocs) as well as in the South of the Amazon ba- 
sin. Following the spatial distribution observed in GiEMS 
estimates (Figure 1 ), floodplains associated with major river 
channels are well delineated. 

[ 33 ] Given the absence of other independent, large-scale, 
multiyear surface water volume estimates over the Amazon 
basin, the seasonal and interannual variability of our results 
over 1993 2007 are evaluated by comparison (Figure 9) 
with related hydrological variables, i.e., the in situ river 
discharge measured at Obidos (see Figure 1 for its location), 
basin-scale estimates of precipitation from GPCF [Adler 
et al ., 2003], and the variations of GRACE TWS (for the 
period 2003 2007). 

[ 34 ] Showing a strong seasonal and mtcrannual variability, 
Figure 9a shows that the time scries of the Amazon River 
discharge is closely linked to the total amount of surface 
water volume in the whole basin. The maximum lagged cor- 
relation is of 0.90 (1 80 months is used to calculate the linear 
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Figure 10. Major droughts in !hc Amazon basin during 1993 2007 as scon by the satellite-derived sur- 
face water volume anomalies and compared to in situ Amazon River observations at O bid os. (a) During 
October November 1997 (averaged and relative to the mean 1993 2007). (b) During September 
October 2005 (averaged and relative to the mean 1993 2007). (c d) Corresponding 3 year anomaly time 
series of Amazon bas in-scale surface water volume (black) and in situ river discharge at O bid os (red), 
(c f) Corresponding 3 year anomaly time scries of surface water volume (multiply by 1 0 for visual purpose) 
aggregated over a large region [0°S 4°S; 57° W 60° W] west of Obidos location (black) and in situ river 
discharge at Obidos (red). 


correlation coefficient, giving p < 0.01, with R > 0. 14) with 
tire surface water volume preceding by l month the in situ 
discharge. On average (for 2003-2007), as shown by the 
climatology in Figure 9b, the increase and the maximum in 
surface water volume often precede the river discharge by 
2 months. Obidos gauging station is located cast of the basin, 
and parts of the water volume variations arc first due to 
variations in upstream locations in the western reaches of 
the basin. The extensive floodplains first store and then re- 
lease large amounts of water to the main rivers and conse- 
quently delay flood waves observed several days/months 
later in Obidos. 

[ 35 ] For the entire Amazon basin and for the 1993 2007 
period, the lagged correlation between the surface water vol- 
ume and the basin-averaged precipitation reaches 0.86, with 
rain preceding the surface water volume by 2 months, and 
/£ = 0.58 with one month lag for the dcseasonalized anomaly 
(computed by subtracting the 1 5 year mean monthly value 
from the monthly time series, Figure 9c). The temporal pat- 
terns between the three variables in Figure 9e arc in good 
agreement and follow alternatively wet and dry events asso- 
ciated with the El Nino/La Nina phenomena. The good agree- 
ment in terms of seasonal and intcrannual variability between 
all variables is confirmed in Figure 9d during the 5 years 


(2003—2007) of overlapping data. The in situ river discharge, 
the surface water volume around Obidos (purple line), and 
GRACE TWS reveal the large anomalous event of the 
2005 drought. 

[36] The droughts that affected large areas of the Amazon 
basin in recent years arc amongst the most severe ones in 
the past hundred years {Marengo et aL, 2008a, 2008b; 
Chen et al., 2009; Tomasella et al. y 2011; Coelho et al ., 
2012; Frappart et al., 2012, 2013] with the 1997-1998, 
2005, and 2010 events still considered as the most excep- 
tional ones in the last 40 years, affecting extensive areas of 
more than 2 million km 2 . The impact on the Amazon rainforest 
during these events was strong, highlighting its vulnerability 
to extreme drought conditions [Philips et al., 2009; Sevan 
et al ., 2009]. Such events have also large potential impacts 
on regional biogeoehemieal and carbon cycles [Philips et al., 
2009; Potter et al., 201 1; Davidson et al., 2012], as well as 
important consequences on human activities and economy. 

[ j?] For the fust time, the signatures of such phenomena on 
the dynamics of surface freshwater volume can be quantified 
directly from observations and analyzed using our new data 
set (Figure 9 and Figure 10). Time series in Figure 9a show 
the associated negative anomalies, with a large deficit in 
terms of surface water volume and river discharge for the 


11,961 


PAPA ET AL.: AMAZON SURFACE FRESHWATER STORAGE 


years 1997 and 2005. The spatial and temporal patterns of 
the droughts during these yearn arc further illustrated in 
Figure 10. The drought of 1997 (Figure 10a, here October- 
November, the drought being in absolute value at its maxi- 
mum during these two months) is characterized by large 
regions of negative anomaly in the south (Mamorc), North 
(Negro and Japura), and East (Amazon) of the basin in good 
agreement with observations of river discharges [Tomasella 
et al. , 2011]. The severe regional drought in the southwestern 
Amazon in 2005 (here September October) is also well cap- 
tured (Figure 1 0b) ? in good agreement with the spatial extent 
as observed in the GRACE data [Frappart et al ., 2012]. For 
both years, the monthly evolution of surface water volume, at 
basin scale (Figures 1 0c and 1 Od) and for the region around 
the gauging station Gbidos (Figures 10c and lOf), show a 
close correspondence with the in situ river discharge, indicat- 
ing a similar deficit for the same period. During October 
November 1997, the amount of surface water stored in the 
entire Amazon basin was about -230 km 3 (or 38%) below 
the October November 1993 2007 average. A deficit of 
-210 km 3 (or 48%) for September October 2005 is found, 
a little bit less than half the GRACE-observed TWS defi- 
ciency of -500 km 3 as given by Chen et al. [2009]. In 
the region west of Obidos [0°S 4°S; 57°W 60°W], these 
deficit values can reach up to 60% in 1997 and 73% in 
2005, in good agreement with the estimates (70%) from 
Frappart et al. [2012]. Investigating the causes and impacts 
of these anomalous events in the Amazon is far beyond the 
scope of this paper, but the new availability of these long- 
term continuous estimates of surface water volume will help 
such future studies. 

5* Conclusion and Perspectives 

[lx] This study presents the first observation- based data set 
that quantifies for the entire Amazon basin the monthly distri- 
bution and variation of surface freshwater storage at -25 km 
sampling intervals over more than a decade, 1993-2007. 
The method is based on a hypsographie curve approach com- 
bining surface water extent from a multi-satellite technique 
(GJEMS) with topographic data from ASTER GDEM. It is 
a three-step process where, to summarize, cumulative distri- 
bution functions of elevation and surface- volume profiles 
arc derived from ASTER GDEM and combined with the 
monthly variations of pixel water coverage from GlEMS in 
order to derive surface freshwater volume. Prior to the last 
step, an evaluation of ASTER GDEM hypsographie curves 
is performed with the hypsographie curves derived from 
EN VlSAT-radar altimeter. And an adequate correction is ap- 
plied to ASTER GDEM database to avoid water volume 
overestimations. The new estimates of monthly surface fresh- 
water volume show realistic spatial structures. The temporal 
variations that arc associated show a strong seasonal cycle 
and interannua! variability. The basin-scale mean annual 
amplitude of -1 200 km 3 contributes to about half of the var- 
iations of GRACE-derived TWS. Monthly surface water 
storage variations for 1993 2007 are evaluated against other 
related hydrological variables such as in situ river discharge, 
precipitation and GRACE observations, showing that the 
seasonal and interannual variations agree well within all vari- 
ables. Finally, we show that the new data set captures well the 
two extreme droughts observed in 1997 and 2005 and helps 


to quantify the water deficit during that period. The amount 
of surface water stored in the entire Amazon basin during 
Octobcr-Novcmbcr 1997 was about -230 km 3 (or 38%) be- 
low the Octobcr-Novcmbcr 1993-2007 average while it 
reached a deficit of -210 km 3 (or 48%) for Scptcmbcr- 
Octobcr 2005, about half the GRACE-observed deficiency 
of —500 km 3 of TWS in the entire basin for that period. 

[ 39 ] ASTER GDEM, and global satellite-derived OEMs in 
general, such as SRTM30 DEM for instance, show a scries 
of characteristics, artifacts, and anomalies that can cause 
significant problems or errors when used for hydrological 
applications [Valeriano et al., 2006; Yamazaki et al., 2012]. 
Some typical problems include the influence of vegetation 
cover and man-made constructions. In ASTER GDEM, the 
acquisition of elevation data using near infrared spectral band 
can introduce errors due to cloud cover, such as very low but 
dense boundary layer clouds in tropical regions. In areas 
of low relief, these small deviations from the true surface 
elevation can cause errors that are difficult to filter when 
representing the true river and floodplain profiles. These 
effects may create artificial islands in the floodplain likely 
to introduce inaccurate elevation in the DEM with conse- 
quences on the hypsographie curve technique. Moreover, in 
order to capture correctly the extreme low storage values 
during droughts, ASTER GDEM should have produced 
credible elevation data for those periods at the low end of 
the hypsographie curves. Such information is difficult to 
verify. In order to improve the performance of DEMs for hy- 
drological applications, a few attempts have been proposed 
in the literature to reduce their uncertainty [Wilson et al., 
2007; Coe et al ., 2008]. However, DEM errors remain as 
one of the main sources of uncertainty in understanding and 
quantifying the interactions between rivers and floodplains. 
On a large to global scale, due to the regional variations in 
the type of errors, no method still exists that addresses all 
possible problems. These issues should be investigated in 
future studies. 

[ 40 ] This unique data set of Amazon surface freshwater 
volume variations over 1 5 years (to be updated until present 
in a near future when global land surface cmissivitics become 
available) removes a crucial obstacle to progress on several 
fundamental scientific questions and can now be used to 
improve our understanding of hydrological and climate 
processes in the Amazon region. 

[ 41 ] For instance, as proposed in Frappart et al. [201 la] 
over the Rio Negro, the surface freshwater storage, together 
with soil moisture products, can be used to separate the inte- 
gral GRACE signal into the contributions of individual stor- 
age components and isolate the variations of groundwater 
storage. Such results would also help to estimate for the first 
time the groundwater recharge in the entire Amazon basin. 
Bringing together these variables with river discharge mea- 
surements, rainfall and surface evaporation estimates in an 
integrated approach will also yield improved overall knowl- 
edge of the different components of the Amazon terrestrial 
water budget. For instance, Azarderakhsh et al. [201 1] dem- 
onstrate the importance to integrate the behavior of surface 
water dynamics when investigating and characterizing the 
closure of the Amazon water budget and its link to natural 
climate variability. In addition, this new 1 5 year data set of 
surface water volume over the Amazon represents an unprec- 
edented source of information for future hydrological or 
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climate modeling of the region. In particular, it will help to 
better characterize and understand the processes during 
extreme events such as exceptional droughts and floods. 
Moreover, it will play a key role in the definition and valida- 
tion of future hydrology-oriented satellite missions such as 
SWOT (Surface Water and Ocean Topography). 

[42] Finally, as GIEMS and ASTER arc available globally, 
this study is also a first step toward the development of such 
database at the global scale. A global surface freshwater vol- 
ume data set is crucial to understand the role of continental 
water in the global water cycle, as well as in the present sea 
level rise. Although progress has been made on quantifying 
the two primary contributors to sea level rise, namely thermal 
expansion due to ocean warming and melting glaciers and 
ice sheets, large uncertainties remain regarding the effect of 
changes in continental water storage, despite recent impor- 
tant results from GRACE. The future estimates of global 
surface freshwater volume will give the opportunity to better 
quantify fluctuations in freshwater flux to the ocean, and will 
help better understanding the role of terrestrial water in the 
present sea level rise and variability [Ra milieu et aL , 2008 ; 
Milly et al. , 2010 ; Pokhrel et al 2012 ]. 
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